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1.0 INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

NASA is currently conducting materials research directed 
toward the production of high quality tri-metal solid solution 
single crystals such as mercury cadmium telluride CHg^.^^Cd^Te) 
for use as infrared detectors. The program includes both ex- 
perimentation in space (Spacelab Experiment No. MPS77F069) and 
ground based research to determine growth characteristics and 
configurations suitable for thio class of materials. One of 
the primary growth requirements of solid solutions as compared 
to elemental and compound semiconductors is the necessity of 
extremely high thermal gradients at the growth interface to 
avoid constitutional supercooling. Growth configurations of 
the vertical Bridgman-Stockbarger type appear best suited to 
produce such high thermal gradients. 

This report presents the results of a study in which a de- 
tailed thermal analysis of a cylindrical HgCdTe sample in a 
Bridgman-Stockbarger crystal growth configuration was conducted. 
The study concentrated on the thermal profile, interface shape 
and position, and the thermal gradients at the liquid-solid inter- 
face. This effort's analysis was based on alloys of HgTe and 
CdTe with compositions approximating 20% CdTe, 80% HgTe. This 
composition results in a bandgap suited for the detection of 
10. 6y CO 2 radiation. 

As part of the initial steps toward successful space growth 
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of HgCdTe, this analysis addressed the sensitivity of the sample 
thermal characteristics to important growth parameters, such as 
thermal diffusivities , thermal conductivities, furnace temperature 
profile, ampoule dimensions, and growth velocity. Physical properties 
of the different constitutents as well as thermal characteristics 
are delineated in Appendix A. Alternate growth configurations 
were defined and evaluated to establish the Bridgman -Stockbarger 
system which best promotes optimum growth conditions. 

Of primary interest in this study is the determination of 
the minimal thermal gradient over the growth interface as this 
determines stability in terms of constitutional supercooling. 
Throughout this report, this thermal gradient, located on the 
longitudinal axis of the sample at the growth interface in the 
direction of increasing temperature the liquid) , will be 

refered to as the "thermal gradient" unless otherwise specified. 

The study was carried out for HgCdTe contained in a fused 
silica crucible as this, to date, represents the only acceptable 
material from a chemical and structural standpoint. 

Included in this study effort was the development of the 
numerical techniques and associated computational models necessary 
to analyze the heat transfer process within the sample and the 
Bridgman-Stockbarger boundary conditions. This thermal analysis 
model has been programmed in Fortran V, and is currently 
operational on the MSFC Univac 1100 system. Detailed descriptions 
of the thermal equations used in this analysis are presented in 
Appendix B. 


1.2 STUDY APPROACH 


As a first step to bound the thermal gradient achievable 
and to provide a test for the numerical computations, the so 
called "constant surface temperature" case was analyzed. This 
case considers an infinitely long cylinder of uniform thermal 
conductivity with upper half and lower half surface temperatures 
at fixed values of Tj^ and T^ respectively. Under these conditions, 
the problem can be solved analytically and yields under suitable 
restrictions the maximum thermal gradient achievable at the 
center of the sample and at the growth interface. 

The next step was to consider the realistic configuration 
with the sample enclosed by a quartz crucible and various thermal 
couplings to a furnace wall. Parametric studies were carried 
out to determine the effects of the following: 

e Crucible characteristics 

• Thermal coupling to the furnace wall 

• Thermal profiling in the hot end 

• Cold end temperature 

e Sample radius 

e Crucible thickness 

e Thermal conductivities 

e Growth velocity 

Specifically "the thermal gradient", interface shape, and 
interface thermal stability were studied as a function of the above 
variables leading to specific growth configurations which can in 
principle achieve the thermal characteristics required of the 
HgCdTe system. 







This study led to the following three distinct growth con- 
figurations : 

• Vacuum Coupling: The volume coupling the quartz cru- 

cible to the furnace wall is evacuated. 

• Gas Coupling: The volume coupling the quartz crucible 

to the furnace wall is filled with a gas with conduc- 
tivities ranging from 0.1 to 0.5 watt/meter-degree. This 
range encompasses the thermal conductivities for argon, 
nitrogen, and helium -- three suitable coupling agents. 

e Liquid Metal Coupling: The volume coupling the quartz 

crucible to the furnace wall is filled with liquid 
Gallium on the cold end and either a vacuum or gas on 
the hot end. 

The range of thermal gradients achievable with these con- 
figurations was established together with the interface shape 
and stability. 

1.5 SUMMARY OF RESULTS 

As indicated, the ’’constant surface temperature” case was 
analyzed with a detailed account given in a later section. The 
thermal gradient at the growth interface is given by 

where R is the sample radius. For the constant conductivity material 
in quartz with an interface temperature of "00°C, we have Tj^*1000°C 
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(yield point of quartz) and T^"400^C (to obtain a 700°C inter- 
face temperature). Thus for a sample diameter of one centimeter, 
the gradient is 

Q . 2 jl000^4^j. goo°C/cm 


It should be noted that this result is independent of the thermal 
conductivity of the sample so long as it is uniform. With the 
restrictions that the quartz does not exceed 1000°C and that the 
700°C isotherm be flat, and not considering effects of different 
thermal conductivities of the solid and melt, this value represents 
the maximum achievable thermal gradient for a 1 cm diameter sample. 

In all configurativins studied, the effects of growth velocity 
for values of interest to the HgCdTe system were found to be in- 
consequential. In fact, as shown in detail later, growth velo- 
cities as high as 5 cm/hr produce no significant effect on the 
thermal gradient and only a 1% variation was found for growth 
rates as high as 10 cm/hr. In this summary, we do not further 
consider growth velocity as an instrumental parameter in the 
growth configuration. 


All results present in this summary correspond to a hot end 
furnace wall temperature of 1000°C. 

For the vacuum configuration, the following gradients were 


found for various 

R(mm) 

6 

5 

4 

3 

2 


ample radii: 

Minimum Probable 
G(°C/cm) 

263 

300 

34r^ 

422 

533 


Maximum Probable 
G(°C/cm) 

315 

450 

496 

548 

670 
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The lower values represent the most unfavorable cases of thermal 
conductivities included by the uncertainties in present knowledge 
Similarly, the larger G values represent the most favorable cases 
It should be noted that an increase of up to 33% in gradient can 
be achieved by profiling the hot end temperature in the vacuum 
configuration. This effect, as expected, is not as pronounced in 
the other configurations. 

In the gas configuration, the thermal coupling to the fur* 
nace wall is increased with a corresponding increase in thermal 
gradient. The following typical data is for nominal values of 
the conductivities and for a gas conductivity of k -0.5 watt/ 
meter-degree . 

RCmm) G(°C/cm) 

5 507 

4 643 

2 1280 

For lower gas conductivities, a substantial increase in G is 
realized over the vacuum case. For example, for R"'5 mm and 

k *0.1 watt/meter-degree, a gradient of 443°C/cm is realized. 

S 

The Ga configuration was studied both with vacuum and gas 
coupling on the hot end. The thermal gradients achieved, again 
for nominal values of the conductivities, for the Ga-vacuum case 
are as follows: 

R(mm) G(°C/cm) 

5 404 

4 484 

2 670 
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We note that no significant increase is found, in fact, the 
gradients are not as high as the k *0.5 gas configuration. In- 

o 

creasing the thermal coupling on the hot end by the addition of 
gas results in the following for R*5 mm: 

k„ (watts/meter- degree) G(®C/cm) 

0 404 

.5 647 

1.0 705 

y 

In the Ga-vacuum case, the heat input from, the 1000°C furnace wall 
is limiting the gradient achievable; an increase in thermal 
conductivity in the hot end allows the gradient to increase. For 
higher gradients, a method of providing higher hot end temperatures 
must be found. A growth configuration capable of accomplishing 
this is discussed in the next section. 

1.4 STUDY RECOMMENDATIONS 

From the above results, it appears possible to get HgCdTe 
crystal growth within the Bridgman-Stockbarger configuration. We 
have already noticed the significant increase in thermal gradient 
realized by incorporating a liquid metal solution in the extractor 
end to increase cool -end conduction and drive the sample isotherms 
up into the melt. If we now can alter the configuration on the 
hot end to drive these same isotherms back toward the cooler, ex- 
tremely high thermal gradients will result in the vicinity of the 
growth interface. 
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This configuration includes a fused silica screen (or "sleeve") 
which has been inserted into the hot end of the furnace cavity to 
surround the HgCdTe sample and crucible. Because of the quartz 
transparency phenomenon (see Appendix B ) , (quartz is transparent 
to shorter wavelength opaque to longer wavelengths) the sleeve 
effectively screens out the longer wavelength radiaCion, such 
that the outer crucible surface is no longer heated by radiation. 
This implies that nearly all radiation from the furnace ';all passes 
directly to the HgCdTe ingot. Additionally, a gas flow is intro- 
duced between the sleeve and the ampoule for a "localized" cooling 
effect. Finally, a very hot (e.g., 1400°C) wall temperature 
source is used to get copious short wavelength radiation. 

In this configuration, the melt temperature rises substantially 
above the crucible strain point temperature at the inner crucible 
surface. However, the heat flow in the crucible wall is now 
radially outward -- from the radiation captured by the melt ^ 
the cooling gas. Because of this cooled outer surface, the ex- 
ternal crucible temperatures are below the strain point tem- 
perature, giving mechanical strength to the ampoule. The sleeve 
is protected in much the same way by the gas along its inner sur- 
face. Because of the "greenhouse" effect imposed by the silica 
sleeve, temperatures will exist in the melt that are higher than 
the 1000°C rupture limit of the crucible. 

Results from this study indicate that this greenhouse hot-end/ 
gallium bath cool -end configuration will produce temperatures on 
either side of the interface that provide the required gradients 
to avoid interface breakdown within the HgCdTe solution. In 
addition, minor experimental adjustments (e.g., decreasing the 
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sample radius) might be made to the configuration to further im- 
prove the thermal characteristics within the sample. 

Other configurations appear feasible to some degree. For 
small sample radii and small growth velocities, stable growth 
conditions cam be produced in both a vacuum and gas -filled fur- 
nace cavity. It should be realized, however, that in actual fur- 
nace configurations, there will be some reductions in the thermal 
gradients computed in this study. With this in mind, a maximum 
sample radius of approximately 3 mm should be considered, if a 
reasonable growth rate and transient growth length are to be 
achieved. Vacuxim configurations are clearly less favorable than 
those including a gas -filled cavity. 

From the data presented in this study, it is clear that the 
system physical properties (conductivities, emissivities , etc.) 
are determining factors in the feasibility of a given configuration. 
Once these values are known with better precision, the computations 
performed in this study should be refined. 




2.0 DEFINITION OF THE BRIDGMAN- STOCKS ARGER 
GROWTH CONFIGURATION 

The materials processing configuration addressed by this 
study is the growth of single crystals from the melt in a Bridgman- 
Stockbarger configuration. In vertical Bridgman-Stockbarger 
growth, the sample material is contained in a suitable crucible, 
and this ampoule (sample and crucible) is placed in a furnace 
which provides the desired profile. Figure 2-1 portrays the 
system under investigation. The sample, crucible, and furnace 
are designed to obtain the highest possible degree of cylindrical 
symmetry in order to effect an appropriate solid-liquid interface 
shape. In some cases, the sample is rotated to ensure better 
symmetry. Directional solidification results from: (a) the 

movement of the ampoule through the furnace cavity, (b) the move- 
ment of the furnace surrounding the fixed ampoule, or (c) changing 
the temperature profile with respect to time. The first two methods 
are equivalent for an ideal non-rotating system, while the third 
is equivalent to the first two only with certain thermal profiles, 

t 

over a restricted segment of the sample, and, of course, with 
the proper time dependence. 

In the following paragraphs, the Bridgman-Stockbarger crystal 
growth configuration under analysis in this study is described 
in more detail. This effort's analysis was based on the sample 
material consisting of solutions of the compounds Mercury Telluride 
(HgTe) and Cadmium Telluride (CdTe) . Descriptions of the physical 
properties of these compounds are given in Appendix A. 
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Figure 2-1 Bridgman-Stockbarger Growth Environment 
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2.1 GROWTH OF HgCdTe SINGLE CRYSTALS 


In analyzing growth from the melt» this study is concerned 

virh many aspects, including the type of material to be grown, 

tj.e processing temperature and the containment of the molten 

material. In the case of HgCdTe, we are dealing with a mixture 

of a semimetal (HgTe) and a semiconductor CCdTe) , with melting 

te iperatures of 668°C and 1092°C, respectively.^ The exact molar 

percentage of the constituents is usually specified by writing 

Hg^_^Cd^Te (0 « X ^1) where, as indicated, the molar percentage 

of CdTe is lOOx. Fused silica presently represents the only 

2 

acceptable crucible material; however, due to the high vapor 
pressure of Hg; this places a severe restriction on the maximum 
allowable temperature. Since this limits the achievable thermal 
.gradients, growth rates in turn must be limited to avoid con- 
stitutional supercooling. 

Taking the above-mentioned aspects into consideration, the 
growth conditions and the furnace characteristics necessary to 
achieve these conditions must be addressed. One of the major 
considerations is the thermal gradient at the growth interface 
necessary avoid constitutional supercooling, which leads to 
interface breakdown. Constitutional supercooling is a form of 
supercooling caused by a change of solute-solvent ratio due to 
he preferred incorporation of one of the components in the solid 
In the case of HgCdTe, this condition is severe, arising from the 
larg.' : eparation of the liquidus and solidus curves. Noting the 
phase diagram of the HgCdTe system (Figure 2-2), we see that 
CdTe will be removed from the solidifying melt more rapidly than 
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the HgTe. This means that the melt at some distance removed 
from the interface will be richer in CdTe, and will therefore 
have a higher liquidus temperature. If the thermal gradient is 
not sufficient to assure a temperature above the liquidus, then 
some form of crystallization will be initiated in the liquid 
away from the interface. This effect is clearly aggravated by 
the high segregation coefficient, and may be suppressed only by 
maintaining a high thermal gradient or a slow growth rate. Rapid 
diffusion of CdTe in the melt helps, as it tends to maintain the 
concentration of that component constant up to the interface. 

2.2 CONFIGURATION CHARACTERISTICS 

Several experimenter-controlled parameters and physical/ 
thermal material properties exist which are critical to determining 
realistic temperature profiles within the growth/melt of the 
Bridgman-Stockbarger configuration. As part of this study, a 
literature search of data relating to the B-S technique and 
fundamental data on the densities, specific heats, thermal 
conductivities, eraissivities , etc. of HgCdTe and fused silica 
was completed and synopsized in the following paragraphs. 

2.2.1 Experimenter -Controlled Parameters 

The ampoule/furnace configuration for the Bridgman-Stockbarger 
technique was illustrated in Figure 2-1.. Theoretically, there is 
an inf initely-thin insulated region in the furnace separating 
the hot and cool ends. Temperatures in the hot end can be 
expected to range from less than 1400°C (ferro-nickel furnace 
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elements) to 2500°C (tungsten elements in vacuum) . Cool end 
temperatures in this study ranged from +200°C down to absolute 
zero (0°K, -273°C). Growth rates (i.e., ampoule or furnace 
movement) of from one-quarter to one centimeter per hour are 
realistic; however, growth velocity was parametrically analyzed 
in this study up to one meter per hour. 

The physical dimensions of this system, also independent 
parameters in this study, are typically established relative to 
the size of the sample. To avoid the possibility of ampoule 
rupture of high temperatures (due to the high mercury vapor 
pressure) , the thickness of the crucible is approximately 60% 
of the sample radius. A nominal sample diameter is 10 millimeters, 
making the fused silica thickness 3 millimeters. For this con- 
figuration, the furnace cavity may be 20 mm in diameter, leaving 
approximately 2 millimeters between the inner furnace wall and the 
outside of the crucible. 

The length of the ampoule is not critical to this analysis. 
The behavior of the thermal profiles and gradients within the 
ampoule were evaluated in the region within five sample diameters 
on either side of the liquid-solid interface. 

As was mentioned previously, fused silica presently re- 
presents the only acceptable crucible material for this crystal 
growth configuration at the high temperatures contemplated. 

However, silica too has its limitations, and must be restricted 

o 

to less than 1000 C under the high vapor pressure of mercury, 
else rupture may occur. The optical transmission (i.e., radiation) 
characteristics of fused silica in the ultraviolet, visible 
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and infrared spectra are discussed in detail in Appendix B of 
this report. 

2.2.2 Parameters Representing Physical Properties 

Many of the parameters addressed in this thermal analysis 
represent physical properties of materials > and as such are not 
under the control of the experimenter. In some cases, these 
properties are not well known, or are somewhat variable from 
sample to sample Ce-g*> the thermal conductivity of fused silica). 

As such, it was useful in this analysis to treat these properties 
as if they were controllable to observe the effect of variation 
or uncertainty in these elements in the output thermal profile. 

Such an approach also extends the usefulness of the program to 
systems containing materials other than HgCdTe. 

The physical properties of primary interest in the Bridgman- 
Stockbarger system describe the thermal characteristics of HgCdTe 
(the sample) , fused silica (the crucible) , and the coupling agent 
between the crucible outer surface and the furnace wall. These 
thermal properties include the thermal conductivities, emissivities , 
and diffusivities for the sample, crucible and coupling agent 
materials. A detailed discussion of these material properties 
is included for reference in Appendix A of this report; study 
assumptions regarding these properties are also presented. 

2.3 GROWTH ALTERNATIVES 

From the available configuration data, several interesting 
growth configurations have been established for detailed analysis 
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under this study. These configurations define different furnace, 
cooler, crucible and cavity arrangements for the crystal growth 
process. Inherent in these different configurations are a 
wide range of requirements relating to laboratory equipment, 
materials and hence cost. 

The assumed baseline configuration is a uniform temperature 
heater with a uniform temperature heat eixtractor, with radiative 
coupling. The hot -end input power is limited to that which may 
produce quartz rupture temperatures of 1000°C or greater. Room 
temperature and refrigerated and heat extractors (down to 0°K) 
were considered for the cool end of the furnace cavity. The 
sensitivity of ingot thermal profiles and gradients to this 
cool -end temperature was investigated for temperatures ranging 
from -273°C to +200°C. 

Alternate heater end configurations were investigated in an 
attempt to optimize the thermal gradient at the liquid-solid 
interface. The maximum allowable quartz temperature limitation 
was enforced, but various temperature boosts near the interface 
on the heater end were specified, and the payoffs evaluated. 

By specifying several maximum temperatures at the interface with 
decaying thermal profiles along the heater wall such that at no 
point the outer crucible surface exceeded 1000°C, numerous state- 
of-the-art furnace configurations were modeled: 

• Ferro-nickel (e.g., .NICHROME, KANTHAL) heater elements -- 
providing up to 1400°C temperatures close to the interface. 
Possibly the least inexpensive furnace configuration for 
crystal growth. 
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• Platinum alloy (e.g., Pt-401 Rh) heater elements -- 
another inexpensive furnace design in which temperatures 
reaching 1800°C in air can be achieved with the platinum 
windings . 

• Tungsten heater elements -- capable of achieving 2500°C 
in a vacuum; here, however, a quantum jump in cost is 
incurred because of the machinery associated with creating 
the vacuum. 

Thermal analyses of these configurations were performed, 
and the payoffs in thermal gradients analyzed with system cost 
effectiveness in mind. If the relative increase in thermal 
gradient, with a consequent increase in crystal growth rate, 
does not justify the corresponding cost or complexity of the 
system configuration, then perhaps multiple systems producing 
an overall higher rate is the rational approach. 

Other furnace-crucible arrangements were considered. The 
effect of sample radius (from 2 to 6 millimeters) within the 
same size furnace cavity was assessed. The crucible thickness- 
tb-saraple radius ratio was checked at 60% for all radii and, 
in addition, was reduced to as low as 33% for some samples. 

Finally, for various hot and cold end temperature profiles, 
the furnace cavity was assumed to be filled with various gases 
(e.g., argon, dry nitrogen, helium), as well as evacuated. 

Thermal conductivities of the cavity medium were varied from 
0 - 0.5 watt/meter-degree. 
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3.0 BRIDGMAN- STOCKS ARGER THERMAL ANALYSIS 

As the objectives of this study indicate, emphasis 
in this analysis was placed on investigating the behavior of the 
thermal profiles and gradients within an HgCdTe sample for selected 
Bridgman-Stockbarger growth configurations. Experimenter-controlled 
parameters (e.g., crucible dimensions, input temperature profiles) 
were analyzed parametrically to establish optimum growth con- 
figurations. In the following sections, the effects of the following 
variations in the B-S growth configuration will be presented: 

• Alternate inner furnace wall temperature profiles --at 
both the heater and extractor ends. Constant surface 
temperatures and boosted (and "deboosted") profiles were 
analyzed. 

• Radiation cooling - vs - liquid bath cooling configurations. 

e Ampoule (crucible/ingot) dimensions. 

e Thermal conductivities of crucible, ingot (melt and solid), 
and ampoule surroundings. 

c Growth velocity. 

As part of this analysis, certain baseline configurations have been 
established which are used as a point of comparison to illustrate 
the effects of the parametric variations being considered. These 
baseline configurations are appropriately described in this section 
also. 

Several special cases of the vertical B-S technique also 
were of interest to this study. These cases correspond to simpli- 
fications in the boundary conditions such that the system differential 
equation presented in Appendix B can be integrated analytically. 
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The analysis of one such case -- The Constant Surface Temperature 
Case --is initially presented in this section to provide insight 
into the physical problem and to illustrate how the thermal pro- 
file is influenced by the physical parameters at the discretion of 
the experimenter. Results from this special case analy is were 
also used to corroborate the numerical techniques used to solve 
the non-integrable problems in this study. 

3.1 CONSTANT SURFACE TEMPERATURE SPECIAL CASE ANALYSIS 

To put the analysis of the Bridgman-Stockbarger growth con- 
figuration into proper perspective, simplifications in the system 
boundary conditions can be introduced to establish limits, trends, 
and deficiencies of the various growth configurations. The 
analytic solution of the Constant Surface Temperature Case is an 
example of the kind of solution which can be used. 

Consider, as in Figure 3-1, the boundary condition 


T • T^ at r ■ R for : > 0 
T « Tj^ at r • R for 2 0, 
where R is the sample radius. 


Figure 3-1. Prescribed Sample Surface 
Temperature Boundary Condition 
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This case approximates the conditions which produce a maximum ther- 
mal gradient at the solid-liquid interface with a maximum allowable 
temperature. Several aspects of this case are interesting. This 
problem for zero growth rate reduces to a solution of Laplace's 
equation with prescribed boundary temperatures. We first note 
that the thermal profiles are. not dependent upon the thermal proper- 
ties, provided the conductivity is given the same value at both ends 
of the sample. The solution is 



Jo 



z>0 


^ ■ 1 



z: 

n*l 


■oJl(R-n) 


e 


4>c 


n 


2<0 


where are the positive roots of (R«) ■ 0. While these 
series provide a formal solution, certain problems arise. In 
the HgCdTe system, one of the most significant results is the 
minimum thermal gradient at the solid-liquid interface. If we 
attempt to compute this value, we note that the series 
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does not converge at r ■ 0. Therefore, this "standard textbook" 
solution is of little use to us in the region of most interest. 
However, it will provide, in the limit, an indication of the 
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maximum possible thermal gradient on axis at the interface as 
follows. 

From the above normalized solutions, the range of is from 
0 to 1, with the midplane value O.S. The cylindrical coordinates 
are r and z, R is the crystal radius, and is the n^^ zero of 
J^. On the axis, r * 0 and >JoCr«j^) » 1. Figure 3-2 presents 
the axial temperature and thermal gradient for this constant 
surface temperature case with zero growth rate. 

In the computer program that produced these results , the 
coordinates used are in units of R so that: 


z' = z/R, r' = r/R 


and the on-axis function is 


«(r’=0) * ^ 


n*l 




♦ « z 
n 


We expect that at z' = 0,$* 0.5, since continuity requires that 
the two equations for (j must yield the same result; therefore, 
<!>(2'=0) = 1 - (z' = 0). The fact that the computer does not pro- 

duce this result is due to the singular behavior of the series 
near z'*0. For all other values of z' , <^(z'*a) * 1 - 
Figure 3-3 illustrates the behavior of <t> for larger values of 
z', indicating that ♦-►I when the axial distance from the inter- 
face is approximately three times the crystal radius. 

The gradient of 4> is taken by dividing the difference between 
successive values of (p by the increment in z'. Since the values 
of p nearest z’*0 are in error, the gradient of <p is likewise in 
error at that point, and, in fact, its sign appears reversed. 



3-4 




The 





0 


Figure 3 


o.t i.: 1.6 2.0 :.a :.3 

:• 


3 Temperature Profile Trend Away from the Interface 
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This is purely an artifact of the computation and illustrates 
the necessity of a careful examination of the convergence of such 
series representations. It appears that approaches the value 

2/3 at the origin, and that this then represents the limiting value 
of the gradient. This limit cannot be reached in real cases, as 
any real boundary conditions will be less favorable than these 
fictitious idealized ones. 

In practical terms, if Tj^ * 1000°C, limited by the mechanical 
properties of fused silica, and if the solidus temperature is 
700°C (the mid-plane temperature), then T^ = 400°C. Thus, Tj^ - 
T^ = 600°C, and for R = 1/2 cm 



Thus, 800°C divided by the crystal diameter appears to be an 
absolute limit on the achievable minimum gradient. The most 
general expression would be 



where T^^^ is the midplane temperature ar.d d is the diameter of the 
specimen. 


3.2 PRELIMINARY CONFIGURATION ANALYSIS 

Consider again the Bridgman-Stockbarger growth environment 
of Figure 2-1. The initial growth configuration considered by 
this analysis was a constant surface temperature profile for 
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both the inner furnace and cooler walls. The hot end temperature 
was specified to be 1000°C (limited by the high vapor pressures of 
mercury and the mechanical properties of fused silica). The extrac- 
tor wall was specified to be -270°C (approximately absolute zero) 
to maximize the thermal gradient at the liquid/solid interface for 
this configuration. (One must recognize the experiment cost impli- 
cations associated with this cooler configuration, however, when 
assessing the payoffs of this imposed external temperature gradi- 
ent. As this analysis will indicate, heat extractor temperatures 
below '100°C have little influence on the system performance.) 

For this configuration, the following additional specifications 
were selected: 


• Thermal Conductivities 


Fused silica : 
Sample (liquid) : 
Sample (solid) : 


2.5 watts/m°C^ 
1.0 watt/m°C 
1.0 watt/ra°C 


(These conductivity values are later varied parametrically 
to assess their individual effects on growth conditions. 

For this preliminary configuration, equal conductivities for 
the melt and solid were assumed to force the thermal gradient 
immediately on either side of the liquid-solid interface to 
be equal.) 


• Ampoule Dimensions 

Sample ingot : 1 cm diameter 

Crucible thickness: 601 of sample radius, or 3 mm. 


t Emissivities 


Fused silica: 

1.0 

Furnace : 

1.0 

Sample : 

0.3 


1 


See Appendix A 
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In addition, the furnace cavity is evacuated; that is, there is no 
conduction between this surface and the furnace wall. Finally, 
there is no ampoule motion, i.e., the growth velocity is zero. 

Figure 3-4 presents the thermal characteristics within the 

sample for this described configuration. The temperature 
profile along both the outer crucible wall and the sample centerline 
(axis) are presented, as well as the axial thermal gradient. Note 
that the outer crucible temperature falls away from the 1000°C 
steady state value closer to the liquid-solid interface, indicating 
that a temperature boost in the heater end near the interface may 
be possible without jeopardizing crucible strength. This 
temperature boost may drive the maximum achievable gradient at 
the interface higher than the indicated 306°C/cm for this con- 
figuration. 

Figure 3-5 presents the isotherms in the ampoule for this 
configuration. These isotherms remain relatively flat for 
approximately 100°C on either side of the HgCdTe solidus 
temperature of 700°C, promoting good crystal growth conditions 
at the interface for high thermal gradients. 

5.3 TEMPERATURE BOOST GROWTH ALTERNATIVE 

■ As the previous configuration analysis indicated, increased 
thermal gradients may be achieved by boosting the temperatures 
of the inner furnace wall near the interface. This configuration 
will yield an exponential temperature decay from a maximum 
temperature (near the interface) to the steady state temperature 
of 1000°C away from the interface. The form of this temperature 
profile (normalized) is (See Appendix B) ; 
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Constant Sur 
Tj^»1000°C 
T^--270°C 


Distance From Solid/Liquid Interface (millineters) 




♦wall * 


where A and B are constants representing the maximum temperature 
and temperature decay rate, respectively. 

For this analysis, boost temperatures from 1300® to 2400®C 
were considered. This range is representative of furnace con- 
figurations utilizing heating elements of ferro-nickel (1200- 
1500°C) , platinum-alloy (1600°-1800° in air) , and tungsten (up 
to 2500® in vacuiom) composition. For each temperature case, various 
decay rates (i.e., values of B) were investigated to determine 
the "critical value" of B such that the temperature profile on 
the outer crucible surface approached (but nowhere exceeded) 
the constant surface temperature limitations of 1000®C. Figures 
3-6 and 3-7 present the thermal profiles for the minimiun and 
maximum TBOOST temperatures considered -- 1300°C and 2400°C, 
respectively. The trend of outer crucible temperature for 
varying decay rates is evident in either case. Also illustrated 
is the axial thermal gradient for the "critical" decay rate, the 
optimum gradient achievable for each TBOOST configuration under 
the physical limitations imposed by the crucible. 

Table 3.1 presents the maxim;am thermal gradients achieved 
for all temperature boost cases analyzed. Note that up to 30% 
increase in gradient (close to 400®C/cm) is realized over the 
previous constant surface temperature case. However, as the 
table indicates, this payoff is limited; as the maximum temperature 
is further increased, the maximum i,.'adient moves away from the 
interface and begins to decrease. Figure 3-8 illustrates this 
trend as it compares the T ■ 2400®C boost thermal profiles with 
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■■1300 C (Maximum Hot End Temperature) 



Figure 3-6 Thermal Profiles for the Minimum Temperature Boost 
(T„=1300°C) Configuration 







TABLE 3.1 

TEMPERATURE 

BOOST EFFECTS 

Tj t'C) 

A 

^CRITICAL 

MAXIMUM AXIAL GRADIENT, 
4T/42 (OC/cm) 

1300 

0.23566 

0.40 

340 

1600 

0.47133 

0.40 

350 

1800 

0.62844 

0.40 

380 

2000 

0.78355 

0.45 

395 

2400 

1.09976 

0.70 

380 


Note: Tg ■ Maximum Boost Temperature along Inner Furnace Wall 

Tj - (T^ - T^) * »B * Tc 
-Bz 

where (fig ■ 1. ♦ Ae (Normalized) 

2 ■ Distance away from Liquid-Solid Interface 
T ■ Cool End Temperature ■ -270°C 

w 
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the constant surface temperature data presented previously. 

Figure 3-9 illustrates the effect of boosting the inner 
furnace wall temperature near the interface on the ampoule iso- 
therms. Isotherms for the TMAX - 1000° (constant temperature 
case), 1300°, and 2400°C are compared. The movement of iso- 
therms toward the interface and into the solid end is evident. 

In all cases, the isotherms in the 600°-800°C temperature range 
remain relatively flat. 

3.4 TEMPERATURE "DEBOOST" ANALYSIS 

A similar thermal analysis was performed for the cool end 
of the configuration to establish the extractor wall requirements 
as a function of temperature differential (hot end - vs - cold end) 
for the quartz out^jr surface. Figure 3-10 presents lines of con- 
stant minimum extractor temperatures required to support the outer 
crucible temperature difference ^ch-^cc for various T^^'s. 

Configuration temperatures were varied around the approximate 
growth faces expected for HgCdTe (700-800°C) . As an example, 
if the hot end quartz surface temperature is ideally 1000°C 
(constant) , the extractor would be able to withdraw only enough 
heat within one millimeter from the interface to lower the quartz 
temperature to 905°C -- and it would take a 0°K temperature at 
the extractor wall to achieve this reduction. Similarly, if the 
cool-end quartz temperature is 600°C, the best temperature differential 
along the quartz that the extractor can support is 50”C. Relaxing 
the extractor wall requirements to more reasonable (i.e., warmer) 
temperatures drives T^j^ down for a specified or, for a specified 

^ch'Tcc, drives both T^.^. and Tj^^. higher, 
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Lowest Kxtractor Temperatures for Various 

T and (T , -T ) Temperature Differentials. 
CC C tl cc 

T and T . are cold and hot end temperature 
cc ch 

on araroule outei surface 1 mm grid assumed 
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Of particular interest is the c lost insignificant temperature 
differential gained (less than 10°C for almost all T^^.) by lowering 
the minimum extractor temperature from +200®C to 0°K (absolute 
zero) ** a range in temperature which has definite cost impli- 
cations on extractor design. This data implies that the tem- 
perature profiles are insensitive to radiative -cooled extractor 
temperatures. These radiation-cooling effects are analyzed in 
more detail in the next section. 


3.5 R.\DIATI0N COOLING EFFECTS 

To further investigate the effects of the radiative-cooled 
furnace configuration, we considered both the 1000°C (constant) 
and temperature-boosted (TMAX*1300°C) hot-end furnace profiles 
with variable constant extractor temperatures ranging from +200° 
to -200°C. (The absolute zero temperature cases were discussed 
earlier.) 

Figure 3-11 presents the axial temperature and thermal 
gradient for the constant temperature case for the specified range 
of extractor temperatures. As the plot indicates, less than a St 
(-12°C/cm) increase in maximum thermal gradient is achieved by 
dropping the extractor temperature from 200°C to 0°C. Table 3.2 
provides the achievable gradients for both hot -end temperature 
configurations over a wide range of extractor temperatures -- in- 
cluding higher extractor temperatures possibly associated with 
crystal growth in solid solutions of other materials. The table 
indicates virtually no payoff in requiring an extractor or cooler 
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TABLE 3.2 EFFECTIVENESS OF RADIATION COOLING 





to produce temperatures less than room temperature (25®C) , viith 
very little loss for extractor temperatures considerably warmer. 

Figure 3-12 presents several isotherm plots to compare the 
effects. of radiative cooling for the two hot-end configurations. 

The shift in isotherms within each configuration is minor for the 
corresponding extractor temperatures. A larger temperature move- 
ment is introduced by the previously-discussed temperature boost 
for these extractor temperatures. 

3.6 EFFECT OF AMPOULE DLMENSION VARIATIONS 

This analysis addressed the effect of ampoule thickness on 
the thermal profiles within the sample. For this study, the critical 
temperature boost (Tj^ * -12S0°C) was assumed on the hot end, 
and a constant extractor temperature (0°, 200°C) was selected. 

The baseline sample radius of 5 mm was varied between' 2 and 6 mm, 
with the crucible thickness adjusted accordingly to maintain 
a thickness-to-sample radius C^C/R) ratio of 0.6. Additionally, 
an ampoule configuration consisting of a 6-mm sample radius and 
a 2-mm crucible thickness (aC/R ■ 0.333) was also analyzed. 

The axial temperature profiles for the 0°C extractor case 
are shown in Figure 3-15 for the minimum and maximum radii con- 
sidered. For the smaller samples (and thinner crucibles) , the 
temperatures are naturally higher in the hot end and cooler in 
the extractor. Accordingly, a higher thermal gradient at the 
interface exists for the smaller diameter samples. Figure 3-11 
illustrates this trend. The same input temperature profile pro- 
duces a thermal gradient for a 4-mm diameter sample that is 


3-23 





1400 



Figure 3-13 Effects of Sample Radius on the Ampoule Centerline Temperature 





over twice that incurred within a 12 mm diameter sample. Table 
3.3 presents the maximum thermal gradient for radii between 2 and 
6 mm, for extractor temperatures of 0°C and 200°C. An inverse 
relationship between radius and thermal gradient would be con* 
sistent with the approximation developed for the special case 

of Section 3.1, i.e. , 

II « ® 

3z ’ y 

Further, Figure 3-15 illustrates that the 3T/3z - vs - 1/R 
relationship is nearly linear for the range of radii considered. 

As the curve indicates, 3T/3z increases at a less than constant 
rate by decreasing the sample radius, even though the crucible 

wall thickness and furnace diameter are kept proportional to R. 

Figure 3-16 presents a comparison of the isotherms for the 

R*2 nim and R=6 mm cases. The higher axial gradients associated 
with the smaller sample radius are again evident by the ’’compressed” 
isotherms, i.e., a greater change in temperature over the same 
axial distance. Hence, for R*2 mm, the hot-end temperatures 
are hotter and the cool-end temperatures cooler, as previously 
indicated in Figure 3-13.. 

Analysis thus far has assumed a constant crucible thickness - 
to-sample radius ratio (aC/R) of 0.6. If this ratio were de- 
creased within the same overall system geometry Ce^g** a 6-ram 
radius ingot in a crucible 2 ram thick) , an increase in thermal 
gradient is realized (Figure 3-17). This improved gradient 
must be carefully weighed against the possibility of crucible 
rupture at these high temperatures introduced by the thinner' 
container . 




TABLE 3.0 EFFECTS OF CHANGE IN SAMPLE DIAMETER 


TEMPERATURE -BOOSTED HOT END 


FURNACE TEMPERATURE SAMPLE RADIUS 

HOT END(Tg) COOL END (T^) R, ram 


•1250°C 0°C 6 

5 

4 

3 

2 


-1200°C 200°C 6 

5 

4 

3 

2 


MAXIMUM THERMAL 
GRADIENT (4T/a2) 


280°C/cra 

3150c/cm 

385°C/cm 

455°C/cm 

585°C/cm 


270°C/cra 

300°C/cm 

350°C/cm 

425°C/cm 

53S°C/cm 




radient - versus - Sample Radi 









Figure 3-16 compares the isotherms within the two 6 -mm 
radius samples. The higher thermal gradients inherent in the 
thinner crucible case are again evident around the interface 
(Figure 3-16c). 

3.7 THERMAL CONDUCTION EFFECTS 

The thermal conduction characteristics of the various materials 
composing the Bridgman-Stockbarger configuration play an important 
role in the heat transfer process required for crystal growth. 

As discussed previously, relatively little direct measurement 
data for thermal conductivities of HgCdTe alloy systems is available. 
Because of this uncertainty, it is worthwhile to parametrically 
investigate the thermal profiles within the samples as a function 
of the assumed thermal conductivity of both the solid and the 
melt in the HgCdTe system. From the considerations presented 
in Appendix A, we have adopted nominal values for the thermal 
conductivity of HgCdTe of 0.5 watt/meter degree and 0.8 watt/meter 
degree in the melt and solid, respectively. 

The conduction characteristics of fused silica are more pre- 
cisely known (see Appendix A. 2.1); however, conductivity of quartz 
is significantly influenced by impurities in the fused silica, as 
well by its temperature history. To investigate the implications 
of these effects, we have varied the crucible conductivity between 
1.0 and 2.5 watts/meter-degree to assess the sensitivity of the 
sample thermal profiles to this conduction medium. 

The thermal coupling agent between the cruciole and the fur- 
nace wall is extremely instrumental in achieving the required 




thermal gradients for growth of HgCdTe crystals. This is parti- 
cularly true in light of the relative insensitivity of the sample 
thermal condition to radiation cooling and the inability to 
attain high gradients via this method. Various gasses are 
available for use as this coupling agent, having thermal con- 
ductivities varying from 0.06 watt/meter degree (dry nitrogen) 
to 0.36 watt/meter degree (helium). For this analysis, we have 
considered a range of conductivities from 0 (vacuum) to U.5 
watt/meter degree. 

One final configuration was also investigated -- the con- 
sideration of liquid metal (specifically, gallium) as a heat 
transfer agent at the cold end of the system. The introduction 
of this cool-end bath produces a thermal conduction capacity of 
~34 watts/meter degree --a much better heat transport capability 
- than radiative cooling. 

Table 3.4 summarizes the effects of the system thermal con- 
ductivities on the thermal gradient along the sample centerline. 

An increase in this gradient is generally produced by higher 
solid or gas conductivities, and by lov:er melt or quartz con- 
ductivities. The isotherms of Figure 3-18 support these ob- 
servations. For example, a lower melt conductivity inhibits rapid 
heat transport toward the interface, which effectively compresses 
the sample isotherms and increases the gradient in the hot end 
(compare Figure 3-18a and 3-18b). Similarly, if the solid thermal 
conductivity is increased, the solid-end heat gradient decreases 
(see Figure 3-18c), which causes a gradient build-up on the liquid 
side of the interface. Figures 3-18 (d) and (e) indicate the iso- 
therm movements for decreased quartz conductivity and increased 





TABLE 3.4 EFFECTS OF THERMAL CONDUCTION 


• BASELINE 

(k -2.S.1C -0.8, 
kJ-0.5,k|-0) 


• VARYI.NG CRUCIBLE 
CONDUCTIVITY, k^ 


MAXIMUM GRADIENT (G) 


400°C/cm 


EFFECT 




415”C/cm 

42S®C/cm 

445°C/cm 

475°C/cm 

500°C/cm 



Decreasing k results 
in increasing G 


• VARYING SOLID 
CONDUCTIVITY, k. 


• VARYING MELT 
CONDUCTIVITY, k. 


« • VARYING GAS 
COND'. :tivity, k 




540°C/cm 

515°C/cm 

475°C/cm 

400°C/cm 

370°C/cm 




380°C/cm 

425°C/cm 

445°C/cm 

505°C/cin 

580°C/cm 




545“c/cm 

620°C/cm 


715^C/cra 

740°C/cin 
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gas conductivity, respectively. Figure 3-18(e) particularly 
illustrates the compressed isotherms high gradients) 

associated with increased thermal conduction through the coupling 

agent between the crucible and the furnace wall. 

As a means of increasing the thermal gradient, a liquid gal- 
lium bath at one or both ends presents an attractive option. Since 
gallium has a low melting point (29^C) , a low vapor pressure, and 
a good thermal conductivity (34 watt/meter-degree at 60°C) , it is 
compatible for use in an evacuated system and is easy to work with. 
Figure 3-19 presents axial temperature and thermal gradient pro- 
files, and Figure 3-20 the isotherms for the case of a gallium 
bath used for cooling and vacuum at the hot end of the furnace. 

A gradient of 650°C is achieved, though this is not at the desired 
700°C point. Appropriate temperature adjustments could correct 
this with a small loss in gradient. 

As Figure 3-18 (a) indicates, the critical 700° (crystallization) 
isotherm for HgCdTe has been displaced by almost 2 mm from the mid- 
plane position of the B-S configuration for the Tj^*1000°C, T^*0°C 
configuration. To drive the isotherms back toward this mid-plane 
position, the extractor temperatures must be raised accordingly. 
Figure 3-21 illustrates the desired isotherm conditions for both 
the constant surface temperature and the temperature boost con- 
figurations. The required extractor temperatures for these cases 
are 400°C and 330°C, respectively. The corresponding achievable 
gradients are 325°C/cm and 360°C/cm -- indicating an -lOI increase 
in thermal gradient with the boosted temperature configuration. 
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3.8 EFFECT OF GROWTH VELOCITY 


The steady state differential equation for the temperature 

sT 

T within the ampoule includes a heat transport term 6 where 

9 Z 

6 is the growth rate v divided by the thermal diffusivity of the 
material. As indicated in Appendix A, we can express the dif* 
fusivity as 

-7 2 

Diffusivity^ ■ 8.17 x 10 meter /sec 

where is the thermal conductivity of medium i, if we use the 
Dulong and Petit specific heat and a density of 7600 kg/m^. 

For all analyses presented thus far, we have assumed no 
growth rate, i.e.,v*s«0. In this analysis, we now para- 
metrically consider growth rates from one millimeter per hour to 
one meter per hour to investigate the effects of ampoule movement. 
The growth configuration considered was the constant temperature 
furnace wall, with Tj^ • 1000°C, T^ - 0°C. 

Figure 3-22 illustrates the trends of the axial thermal 
gradient for the range of growth velocity considered. Table 3.5 
summarizes the maximum gradients achieved for this range of 
growth velocities. Note that growth rates below -5 cm/hr have 
little, if any, appreciable effect on the maximum gradient near 
the interface. Less than 5% improvement in effectiveness is 
realized with even faster growth rates, until such time that 
the growth rate becomes large enough that the isotherms are 
swept almost completely into the solid end and the gradients be- 
gin to decrease (Figure 3-23). 
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Distance from Solid/Liquid Interface (millimeters) 
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APPENDIX A 

SYSTEM THERMAL PROPERTIES 
A.l PROPERTIES OF HgCdTe 
A. 1.1 Thermal Conductivity 

Very little reliable work hs.s been done in the direct measure- 
ment of the thermal conductivity of HgCdTe alloys , and so we are 
forced to rely on extrapolation and interpolation to a great de- 
gree. We do have some data for CdTe and for HgTe which are help- 
ful . 

4 Q 

Long and Schmit report that at 77 K CdTe has a thermal 
conductivity of 40 watts/meter-degree, and that k varies with 
T'-^'-^ at least to 300°K. Also that HgTe at 77°K has k » 15 watts/ 
meter-degree, anc that k varies with T*^ to 300°K. Furthermore, 
they report "a shallow minimvim" in the alloys for a composition 
of Hgg 2^®’ happens to be close to the composition 

of particular interest. 

Zanio * shows a curve for the thermal conductivity of CdTe 
which indicates that the T*^ law is probably valid to and beyond 
100C°K. 

6 

The TPRC data series gives plots for the thermal conduc- 
tivity of Hg ggCd Q^Te and Hg gjCd gyTe up to 300°K. 

Polycrystalline material i' expected to have a lower thermal 
conductivity than monocrystal, but this effect should not be 
important at high temperature. 

From the above considerations we have adopted single best 
values for the thermal conductivity of HgCdTe in the melt and 
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in the solid of 0.5 watt/meter-degree and 0.8 watt/meter-degree 
respectively, and we have tested values varying above and below. 

We have performed specific analyses using equal values in the 
solid and liquid phases where it was important to show the effects 
of varying other parameters without the interference of a step 
change in conductivity. 

A . 1 . 2 Emissivity 

The emissivity data for HgCdTe is also far from copious, and 
the use of sharply defined values would be dubious in this work 
in any event as much depends on the surface condition of the 
sample and on the nature of the crucible (fused silica) interface. 
It is known, for example, that as solid HgCdTe cools it shrinks 
back from the crucible wall, and that at the instant of separation 
the effective total emissivity drops sharply. 

Based on TPRC ^ data for CdTe for the hemispherical infra- 
red spectral emittance, we have chosen 0.3 as the baseline value 
for the total emittance in both solid and liquid. This value is 
higher than the value typical of materials such as germanium, but 
it is justified. 

A. 1.3 Thermal Diffusivity 

The thermal diffusivity is calculated from the equation 

diff * k/pC 

where k « thermal conductivity 
0 » density 
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C ■ specific heat 


The thermal conductivity is discussed above. 

We have excellent data on the mass density of 
e.g., Blair and Newnham * , who give the data for the following 
table: 


X 

oCg/cm^j 

0 

8.07 

0.05 

7.91 

0.1 

7.84 

0.2 

7.62 

0.4 

7.17 

1.0 

5.81 


It was judged unproductive to make an exhaustive study of 
the specific heat, as the thermal conductivity is so poorly 
known. Thus we have simply used the Dulong and Petit law that 
the specific heat per mole is 3R, We use the classical value 
rather than the Debye law as we are working near the melting 
point, well above the Debye temperature. Hence, C = 25 joules/mole 
degree. For HgQ gCdg 2 ^ 6 . the molecular weight is 0.31058 kg/ 
mole and from the above table the mass density is 7620 kg/meter"^. 
Thus the molar density is 24.6 x lO"^ mole/meter'^ . The diffusivity 
is given by 

Diffusivity ■ 2 X 24.6 x 10^ x 25 * ^ ^ 

where k is in watts/meter-degree . 
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A. 2 PROPERTIES OF FUSED SILICA 


A. 2.1 Thermal Conductivity 

The thermal conductivity of fused silica glass requires a 
special treatment for this analysis > as either of the two pro- 
cesses (radiation photon or phonon) may dominate, depending on 
temperature. Furthermore in a system in which large gradients 
are present, the phonon and photon fields will be very far from 
mutual thermal equilibrium, and care must be exercised in the 
definition of temperature. 

Fortunately, for fused silica, it is possible to make con- 
ductivity measurements which exclude radiation, and photon mean 
free paths are either very short or very long compared to the 
system dimensions, depending on wavelength. Discussion of this 

9 

can be found in Kuigery, Gowen and Uhlman , and the application 
is discussed in Appendix B of this report. 

Data for the spectral photon transparency of fused silica 
have been taken from the Awersil bulletin for T-G commercial 
grade fused silica and approximations were made as xplained in 
Appendix B (Table B.l results). 

Phonon conduction in single phase glasses increases with 
temperature only at low temperatures. At high temperatures, 
the phonon mean free path is limited by the highly disordered 
network structure to the order of interatomic distances. Thus 
above about 800°K the phonon conductivity is nearly constant, 
with slightly higher values for pure silica than for the more 
complex borosilicate or soda-lime g.-..ses. Conductivity is 
significantly influenced by impurities in fused silica, as well 
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as by its thermal history. 

We have taken the fixed value k^ ■ 2.5 watts/meter-degree as 
our phonon conductivity baseline value, and have varied this for 
effect in the parametric analysis. 

A. 3 PROPERTIES OF THE THERMAL COUPLING AGENT BETWEEN THE CRUCIBLE 
AND THE FURNACE WALL 

A. 3.1 Thermal Conductivity 

Three cases were considered. The coupling agent may be a 
vacuum allowing only pure radiative transfer, it may be a gas 
allowing both radiation and conduction (convection is not allowed) 
or it may be a liquid metal such as gallium. 

The various possible gasses have widely varying thermal con- 
ductivities. For example, heliinn at 1000°K has a conductivity 
of 0.36 watt/meter-degree, while dry nitrogen at 1000°K has a 
conductivity of only 0.063 watt/meter-degree. Thus within limits 
the gas conductivity may be considered an experimenter controlled 
parameter. The CRC Handbook gives useful tables for gas con- 
ductivity. 

Liquid metal was considered as a heat transfer agent only at 
the cold end of the system, and for this purpose gallium was chosen 
as an example. The thermal conductivity of gallium at 60°C is 
34 watts/meter-degree, while the melting point is slightly below 
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APPENDIX B 

NUMERICAL SOLUTION TO THE THERMAL PROBLEM 

In formulating the thermal problem for the Bridgman-Stock- 
barger configuration, we did not initially specify the charac- 
teristics of the furnace system. Rather, we specified the thermal 
conditions required within the material sample, together with the 
thermal characteristics and constraints of the B-S system. This 
specified environment defines the boundary conditions and physical 
constants for the problem under analysis. If a unique thermal 
profile is not dictated by these boundary conditions, then it is 
possible to optimize some aspects of the problem -- for example, 
developing a system with minimal power consumption, maximum 
thermal gradient at the interface, etc. Once the desired thermal 
conditions for the ampoule have been determined, the ampoule sur- 
face temperature, the heat flux, and the coupling effects to the 
furnace wall imply the furnace specifications. 

This effort, has been concerned with determining optimum 
ampoule thermal conditions and the generation of furnace speci- 
fics "ions, and did not address the performance of a particular 
furnace design. Thus, the numerical techniques used to solve 
this thermal problem were not concerned with furnace design other 
than the inner wall and the associated thermal coupling to the 
outer crucible wall. 

B.l THE DIFFERENTIAL EQUATION 

In formulating this thermal problem, we have made the 
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following assumptions: 

• There is cylindrical symmetry within the ampoule. 

• All time 'dependent transients have subsided. 

'• The solidification rate is sufficiently slow that the 
latent heat of fusion can be ignored. 


This latter assumption can be relaxed if necessary, but the 
heat of fusion is negligible in most practical systems -- 
especially those requiring very small mean growth rates 
(i.e., solid solutions). 

The steady state differential equation for the temperature 
T within the ampoule then is 


k a‘T ^ 1 /aT\ ^ 3^T 

7 ^ ^ 



( 1 ) 


where the thermal diffusivity K/pC^ is assumed constant for 

any one phase or material (solid, melt, or silica container). 

In this equation, k is the thermal conductivity, p is density, 

C is the heat capacity, r and Z are the cylindrical reference 

coordinates (r is radial position, z is the axial position 

away from the heater-cooler boundary) , and v is the growth 

rate (i.e., v ^ is a heat transport term). 

32 

If we let 


T'T^ 

ttt 


we have normalized the temperature from 0 to 1 as T varies from 
ic to Tj^, the low and high temperature as>Tnptotes of the system. 
'Since d(}i/dx « dT/dx where x is any coordinate, substitution 

-T 

‘h *c 

of for T in equation (1) is immediate so that 
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0 , with 6 ■ 


( 2 ) 



i (o±\ 

r \ar/ 




8 ^ 
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The above differential equation, which reduces to Laplace’s 
Equation when there is zero growth rate for the system, has 
been solved using finite difference techniques to perform 
the Bridgman-Stockbarger thermal analysis. 


B.2 FINITE DIFFERENCE EQUATIONS 

For the finite difference technique, it is necessary to 
describe the physical problem by defining a "grid" or node 
point structure throughout the system of interest, and further 
to specify the interaction of each node point with its neigh- 
boring node points. For the B-S configuration, there are two 
independent variables -- r and z. For this study, we have de- 
fined the region of interest as illustrated in Figure B-1, i.e., 
r > 0 from symmetrical considerations and z extending away from 
the liquid-solid interface in either direction to a distance 
equal to 5 times the sample diameter iz> 0 in the solid, 2 4.0 
in the liquid) . 

This region has been divided into equal divisions such that 
node points are indicated by the coordinates r^ ■ is and Zj • js 
where 5 is the dimension of the node cell. The temperature at 
the ij node then is specified as 41 ^^, that is, the value of ? 
r » r^ and z * Zy To arrive at the finite difference equation 
corresponding to ( 2 ), the following results^® are utilized (ig- 
noring higher order terms) : 
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Figure B-1 Bridgman-Stockbarger Growth Environment 
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d2f 


dx 




m 


For the purposes of this study, it is also desirable to have 


d^f 


dx^ 


^ I M 


m 2m6‘ 




for mj*0 


and 
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Writing these equations for the two variables of interest 
produces 


3 (t 
32 


2 

H 








k <^<fi,j+l ‘ 


IT 


1 Ui,j-i - 2 ^ ^i,j-^i^ 


C3) 

(4) 


+ i it 

r 3 r 

i.j 


- ^ ! (2i-l) - 44j_. , 

2i6^ 1 


C2i^l)^i,l^. 


(5) 


Substitution of these equations into (2) and solving for $ at 
the ij node eventually yields 
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♦l.j ■ 


(1-^) *1*1,3 * 


(1-^) *1,3*1 ^ (‘ * ^)*l,3-l 


i f 0 


( 6 ) 


and 


’O.j 


- 2/3 ^ 






‘OJ-1 


i - 0 (7) 


It should be noted that the effects of the cylindrical system 
have been incorporated into ( 6 ) via the radial correction term 
l/2i. This correction indicates that the temperature at a node 
point is not exactly the mean of the four surrounding node points, 
but is affected more by the larger-radius temperature because 
of the "surrounding" nature introduced by the grid curvature 
effects. This factor becomes less important as the radius 
increases away from the sample centerline. A similar correction 
for a moving ampoule (i.e., non-zero growth rate) is also 
evident in these equations. 

Equations (6) and (7) » representing solutions for internal 
node points, can be used together with the appropriate boundary 
conditions in an iterative fashion to obtain the values of 
(hence T) for all system node points. Care must be taken to 
ensure that the appropriate value of thermal diffusivity corr 
spending to the material at each node point be incorporated into 
these equations. Desired thermal gradients can then be ob- 
tained from the 41 ^^ values within the numerical accuracy of 
the procedure. 
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b.3 treatment of boundary conditions 


Boundary conditions must be treated differently from the 
"internal" node points. Figure B*2 diagrams the Bridgman- 
Stockbarger growth configuration and indicates the types of 
boundary conditions encountered in this thermal problem. 

A fused silica (also called "quartz") ampoule is in a furnace/ 
cooler (also referred to as "extractor") combination. This 
ampoule contains an ingot, the top part of which is melted, 
the lower part recrystallized. Heat flow is into the top part 
of the ampoule, out the bottom. The conditions at the various 
interfaces (e.g., crystal -melt , melt-crucible, crystal-crucible, 
crucible-furnace) are set by laws of continuity. These laws 
are determined by the nature of the particular interface, and 
may generally be approximated by one (or a combination) of 
the following boundary conditions: 


• Prescribed surface temperature or surface temperature 
functionT 

Due to materials limitations (e.g., fused silica rupture 
strength) , it may be meaningful to consider cases where 
the temperature is prescribed on a system surface such 
as the sample, crucible, or furnace wall (Figure B.3). 
While it may not be possible in practice to fix the exact 
temperature of the melt or crystal surface, nevertheless 
the differential equation (2) may be integrated directly 
under these conditions t provide solutions which give 
valuable insight into the nature of the physical problem. 

• Heat conduction across a surface. 


In this case, heat flow continuity requires that 





the Bridgman-Stockbarger System 










where k- and k- are the thermal conductivities on the 
two sides of the surface, and (3T/3n)i represents the 
thermal gradient normal to the surface in region i. 
This case might be appropriate for the molten material 
in contact with the crucible wall (Figure B-4), where 




£ 


with subscript q referring to values in the crucible 
and £ referring to values^in the molten liquid. 


• Radiation -- Liouid/Solid surface in contact vith a trans - 
parent medium. 


Such a boundary condition transfers heat by radiation as 
well as vapor conduction. An example of this boundary 
condition is the coupling between the exterior of a 
crucible and the interior of a furnace (Figure B-5) . 

This radiation flux takes the form S(T. - T- ) , where T^ 
and T2 are " temperatures of the two'^radiating surfaces 
(i.e., the ti.w surfaces coupled by the transparent medium). 
The radiation factor S for the B-S thermal problem includes 
the emissivities of the two surfaces, Boltzman's Constant 
and a transparency fraction for blackbody radiation from 
a temperature source. These factors will be described 
more fully in iollowing discussions relating to the 
radiative boundary conditions specific to the Bridgman- 
Stockbarger growth configuration. 


Using these classes of boundary conditions, the Bridgman- 
Stockbarger growth configuration can be completely defined in 
coniunction with the differential equation (2). Consider, for 
exa.nple, the upper half of the system, i.e., the flow of heat 
to the melt from the heater. The following is happening: 


1) Heat is transferred directly from the heater to the melt 
by radiation from the heater in thos'* wavelengths for 
which the silica is transparent. 

2) Heat is transferred to the silica outer surface from 
the furnace in those wavelengths for which the silica 
is opaque. 

3) Heat is transferred from the heater to the gas (if any) 
by conduction. The gas is considered transparent at all 
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wavele.igths . 

4) Heat is transferred to the silica from the gas by con- 
duction. 

5) Heat is transferred to the melt from the quartz by con- 
duction. (This transfer may go the other way if the gas 
is cooled so that heat is conducted from hot melt to 
cool silica to coir gas.) 

6) Heat is radiated to the furnace by the melt at those 
wavelengths for which the silica is transparent. 

7) Heat is radiated from the ampoule to the furnace at those 
wavelengths for which the silica is opaque. 


Schematically, these transfer processes are depicted in 
Figure B-6. The specified boundary conditions for this B-S 
configuration which relate f> these processes are detailed in 
the following paragraphs. 


B.3.1 Solid-Liquid Interface 


The solid-liquid interface in the B-S configuration is a 
conduction boundary condition in which, for every ij node along 
the interface. 


k 


s 



V 11. 
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For this study, we have assumed no radial heat flow along the 
interface; the impact of this assumption is subject for further 
analysis . 

This boundary condition can easily be expressed in the finite 
difference format such that 4 >ij on the interface can be computed 
from its neighbori.ig i|)ij values. Specifically, for any value 
of i , 
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* 

for j • j at the solid-liquid interface. 


B.3.2 Sample/Crucible Interface (Sample Outer Surface] 


Heat flow to and from this interface is both conductive and 
radiative in nature. Conductive heat transfer occurs between 
the sample (either solid or melt) and the crucible. Radiative 
heat flow from the sample surface, and toward this surface from 
the furnace wall, is also introduced through the crucible be- 
cause of the transparency of the fused silica to shorter wave- 
length (^3.5y) radiation. The boundary conditions for this 
interface then are 



where k is the thermal conductivity, z is the emissivit/ of the 
radiating surface, c is the Stefan-Boltzman Constant, T is the 
temperature of the radiation source, and f^ is the transparency 
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i- 

fraction for black body radiation from a source at temperature 
T^. In these equations, the subscript i represents melt, c ■ 
crucible, h * heater, s • solid, and e * extractor. It should 
be clear that in these boundary condition equations, the tern* 

A 

perature terms T can be expressed in terms of . through the 

A XI 

normalizing relationship 

"^x “ ■ "^0001^ ^i,j "^COOL 




The above equations can also be represented in the finite 
difference format such that 4'^^ at this interface can be solved 
explicitly. However, in addition to ensuring that heat flow in 
the radial direction is governed by these conditions, conser- 
vation of energy dictates that axial conduction effects also 
be included. Thus for the sample/crucible interface we have 
the following conditions: 


Liquid-Crucible : 

■ u.j - -i-i.i 
6 


= k. 






TT- 






i,j ■ ‘^i,j-l^ 


or 


2(kj. k^) 


‘'c *1.1, j * * 




where T^^ » (Tj^q.j. - + '^COOL’ 

and Tj^ = CThot * ^coOL) Wl,J Tcool is speci- 

fied) . 
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Solid-Crucible; 


niCj ♦ k^) - k^ ♦ kj * | k, ♦ Hc pi.jn * Hj-l) 




where and are calculated like and Tj^, respectively 

Both equations (9) and (10) represent quartic equations in (j*. . and, 

^ » J 

as such, can be solved via the Newton-Rhapson technique (or equiva- 
lent) to obtain • at this interface based on the temperatures 

^ » J 

of the neighboring nodes and the inner furnace wall. 

More needs to be said about the transparency function f which 
appears in the radiative terms of these boundary conditions. As 
these equations indicate, heat is transferred directly from the 
heater to the melt in the Bridgman-Stockbarger configuration by 
radiation from the heater in those wavelengths for which the 
silica is transparent. For these same wavelengths, heat is 
radiated to the furnace by the melt. In addition, the furnace 
transfers heat to the silica outer surface in those wavelengths 
for which silica is opaque, and likewise receives radiated heat 
from the ampoule in these longer wavelengths. 

Black body radiation intensity, W, varies with wavelength 
and temperature as indicated in Figure B-7. As the figure 
indicates, the fraction of black body radiated power at wave- 
lengths shorter than x , can be expressed as 


Jo Wdx 

Jo 
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Figure B-8 Fused Silica Optical Transmission Characteristics 
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The American Institute of Physics Handbook chapter on radiometry 

(Table 6k-l, Black-Body Radiation Functions) expresses this fraction 

f for various values of XT, where T is temperature in degrees 

Kelvin. The first radiation integral of the A.I.P. table is 

f ■ 1.3 X 10"® for XT ■ 0.050 cm *^K or 500w °K. From the optical 

trai.smission characteristics of commercial grade fused silica 

presented in Figure B-8, if the ultra-violet cutoff of radiation 

is taken to be 0.2y, virtually total transmission is ensured 

for T • 2500°K -- close to the maximum available output with 

tungsten heaters. At 3600°K and 0.2^, T » 720p °K ■ 0.072 cm °K 

-4 

and f » 3.2 X 10 , less than 0.0321. Therefore, at reasonable 

temperatures there is no significant ultraviolet cutoff of black 
body radiation, and only infrared screening need be considered. 

From Figure B-8, we can assume that in the infrared region, 
quartz has less than 5% first surface reflection loss and re- 
mains perfectly transparent for all wavelengths shorter than 
approximately 3.7)j. Similarly, we will assume that the crucible 
will act perfectly opaque (and black) for x>3.7i. The fraction 
of total black body radiated power which is at wavelengths shorter 
than 3.7y can be obtained as a function of temperature by adapting 
Table 6k-l of the A.T.P. Handbook for fused silica. This is done 
in Table B.l. Fo. * 3.6u and 3.8y, the- appropriate temperatures 
T are presented and correlated with f^, the transparency fraction, 
for the specified value of XT. Figure B-9 graphically presents 
this f^ fraction as a function of temperature for the transparent 
wavelengths . When considering the surface radiation from fused 
silica, the fraction 1 - f can be used as an effective total 

A 

emissivity to a first approximation. 


^ A.I .P. , pp 6-154£f , Table 6k-l 
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TABLE B.l FRACTION OF TOTAL BLACK BODY RADIATED POWER 
AS A FUNCTION OF TEMPERATURE 


XT 

1 

Xq«3 . 6y 
T°C 

Xo»3 . 8w 
T°C 

Transparency 
Fraction, £ (%) 


4.7 

-9.84 

0.032 


60.3 

42.8 

0.215 


115.8 

95.4 

0.782 

0.16 

171.4 

148.1 

1.979 

0.18 

227.0 

200.7 

3.946 


282.5 

253.3 

6.690 


338.1 

305.9 

10.11 

0.24 

393.6 

358.6 

14.05 

0.26 

449.2 

411.2 

18.34 

0.28 

504.7 

463.8 

22.82 

0.30 

560.3 

516.5 

27.36 

0.32 

615.8 

569.1 

31.85 

0.34 

671.4 

621.7 

36.21 

0.36 

727.0 

674.4 

40.40 

0.38 

782.5 

727.0 

44.38 

0.40 

838.1 

779.6 

48.13 

0.42 

893.6 

832.3 

51.64 

0.44 

949.2 

884.9 

54.92 

0.46 

1004.7 

937.5 

57.96 

0.48 

1060.3 

990.2 

60.79 

0.50 

1115.8 

1042.4 

f3.41 

0.52 

1171.4 

1C95.4 

65.83 

0.54 

1227.0 

1148.1 

68.07 

0.56 

1282.5 

1200.7 

70.13 

0.58 

1538.1 

1253.3 

72.04 

0.60 

1393.6 

1305.9 

73.81 

0.62 

1449.2 

1358.6 

75.44 

0.64 

1504.7 

1411.2 

76.94 
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For purposes of modeling this transparency fraction in our 
Bridgraan-Stockbarger thermal analyses, Figure B-9 may be appro- 
ximated by analytically-convenient functions. The first function 
suggested is the straight line 

f^ - CT - 150°C)/1400®C, 

which stays within +3% up to T«107S°C, but must be adjusted at 
both low and higher temperatures. A much improved analytic 
function, which is accurate over a wider range of temperatures 
(particularly at higher T’s) and has no discontinuities in its 
derivative, is 



\ 


f 


X 



( 12 ) 


This curve has its inflection at T ■ a , f = 1/2*'; reasonable 
values which we have assumed for our thermal models are a ■ 
425°C, b - 350°C, and C - 2.5. 


B . 3 . 3 Crucible/Furnace Interface (Outer Crucible Wall) 

Similar boundary conditions exist at the outer crucible wall 
as were specified for the inner crucible surface, namely, 


Crucible - Heater: 


k ^ 

*^c ar . . 


- it 

'g 31' 


1,J 
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Crucible - Extractor ; 


k 


3 ^ 
c 3 r 


i J 


k ^ 
g 3r 


1.1 


Cc.,0 l(l-fe)T; 



where the subscript g represents the gas surrounding the ampoule 
(the remaining parameters are defined as before). Naturally, 
if crystal growth is being accomplished in a vacuum, the con- 
ductive term k. 14- becomes zero. Note that for this boundary 
g 3r 

condition, the energy radiated from the two coupled surfaces 
is that radiation for which the crucible is not transparent 
(i.e., the longer wavelengths). That fraction of energy not 
"passed" through to the ingot ir represented by (1 - f^) , the 
"opaqueness fraction". 

The conservation of heat flow expressed by these boundary 
conditions and the axial conduction along the outer crucible 
surface allows us to specify the following finite difference 
expressions for the crucible/furnace interface: 


Crucible - Heater: 


2(1C, ♦ kj) . nj-i: 




( 13 ) 


r 



I 
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Crucible - Extractor: 




+ e 


j * Vi*iJ 

- (l-fc)T/ I (14) 


where, in both cases, - t*^HOT “ ^C00L^*i,j * ^COOL 


The same transparency function £ (as a function of T ) which 
was described in Section 6.3.2 is utilized in these equations 
to establish the fraction of energy radiated to and from this 
boundary -• indicating the degree of opaqueness of the crucible 
at these temperatures. Again, if the ampoule is in a vacuum 
(k ■ 0.), these equations are simplified somewhat. 

O 


B.3.4 Inner Furnace Wall 

As mentioned previously, the purpose of this study was not 
to address particular furnace designs, but tr determine optimum 
ampoule thermal conditions and establish the associated furnace 
specifications. In this regard, of primary importance is the 
analysis of the inner furnace wall temperature profile and the 
resulting thermal coupling effects with the outer crucible wall. 

To assist in this analysis, an input temperature profile model 
was developed to facilitate parametric evaluations of several 
furnace configurations representing external temperature sources 
to the B-S system. The intent of this model was to provide an 
input temperature profile (normalized) along the inner wall as a 
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function of distance away from the solid/liquid interface, ^ ■ 
f(z), for both the hot and cold ends. The form of this nodel is 


— . 


WALL “ 

1. . 

for 2 < 0 

(hot end) 

(15) 

WALL “ 

C,-D*s 

for 2 > 0 

(cool end) 

(16) 


with the constants A, B, C and D positive numbers. By judiciously 
controlling these input parameters, the inner wall profile can be 
adjusted to represent several furnace configurations 'or example, 
if a constant temperature surface is required at bo: .. the hot 
and cool ends, B ■ U « 0 , and C and A are adjusted to create 
the desired normalized temperatures for all z's. Temperature 
profiles also can be produced which decay exponentially from 
^MAX solid/melt interface to a steady-state temperature 

a specified distance away from the interface. The analysis of 
such profiles determines the thermal gradient payoff achieved 
by using a furnace with a tempe'. itu~e boost capability near the 
heater-cooler interface. Note that ..ce temperature profiles on 
each side of the liquid-solid interface can be controlled in- 
dependently . 

The specified inner wall temperature profilr represents an 
external heat source (and sink) which becomes a boundary 
condition for the Bridgman-Stockbarger configuration. Where 
practical, it is desirable to select the finite difference grid 
system such that this boundary (and all boundaries) :s represented 
by a set of node points within this grid. 
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B.4 ITERATIVE PROCEDURE AND ACCURACY 

Several procedures exist for solving the system of finite 
difference equations with the appropriate boundary conditions. 

For steady state problems, the relaxation process is usually 
..apericr. Util'..ing the correct established relationship (either 
equations (6) and (7) for internal grid points, or the appropriate 
boundary condition equation for a boundary node point (8) - (10) , 
(13) - (16) , a temperature can be systematically calculated for 
every node point ij based on the temperature of the neighboring 
nodes. After a value of 4- • is calculated for each node point 
in the system, successive iterations are performed for the entire 
grid until convergence is achieved. For the physical problem at 
hand, convergence can be formulated by examining the change in 
values on successive iterations , requiring ^SjNEW ~ 

<S. 

The accuracy, stability, and time of convergence of this 
iterative technique in practice depends on the particular boundary 
conditions (i.e., configurations) imposed. This study included 
a parametric variation of to establish the sensitivity and range 
of temperature accuracies, and the corresponding variation in 
number of cycles required for convergence. We have also analyzed 
several techniques to improve convergence including 

• using only values from the preceding iteration for all 
neighboring node temperatures 

• incorporating all new (i.e., from the current iteration) 
node temperatures into the calculations whenever possible 

• providing better initialization values for the grid network. 
For example, initial values of the system thermal profiles 
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are better represented by an arc tangent function ( as a 
function of axial distance from the solid- liquid interface) 
than being constant throughout the hot and cold ends . 
Similarly, when successive parametric cases are being 
executed, the solution to one case provides an excellent 
initialization network for the calculations of the next 
case. 


This study has identified other iterative methods for solving 
the system of finite difference equations presented for the 
Bridgman-Stockbarger configuration, with potential to improve 
the rapidity and accuracy of convergence. We recommend that 
future thermal analysis efforts using the finite difference 
techniques and computer model developed during this study 
address these possible enhancements. 
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